arXiv:1503.02213vl [q-bio.NC] 7 Mar 2015 


Effects of the spike timing-dependent plasticity on the 
synchronisation in a random Hodgkin-Huxley neuronal 

network 


R. R. Borges 1 , F. S. Borges 1 , A. M. Batista 1,2 ’*, E. L. Lameu 1 , R. L. 
Viana 3 , K. C. Iarosz 4 , I. L. Caldas 4 , M. A. F. Sanjuan 5 

1 Pos-Graduagao em Ciencias, Universidade Estadual de Ponta Grossa, 84030-900, 

Ponta Grossa, PR, Brazil. 

2 Departamento de Matematica e Estatistica, Universidade Estadual de Ponta Grossa, 

84030-900, Ponta Grossa, PR, Brazil. 

3 Departamento de Fisica, Universidade Federal do Parana, 81531-990, Curitiba, PR, 

Brazil. 

4 Instituto de Fisica, Universidade de Sao Paulo, 05315-970, Sao Paulo, SP, Brazil. 

5 Departamento de Fisica, Universidad Rey Juan Carlos, Tulipan s/n, 28933 Mostoles, 

Madrid, Spain. 


Abstract 

In this paper, we study the effects of spike timing-dependent plasticity on 
synchronisation in a network of Hodgkin-Huxley neurons. Neuron plastic¬ 
ity is a flexible property of a neuron and its network to change temporarily 
or permanently their biochemical, physiological, and morphological charac¬ 
teristics, in order to adapt to the environment. Regarding the plasticity, 
we consider Hebbian rules, specifically for spike timing-dependent plasticity 
(STDP), and with regard to network, we consider that the connections are 
randomly distributed. We analyse the synchronisation and desynchronisa¬ 
tion according to an input level and probability of connections. Moreover, 
we verify that the transition for synchronisation depends on the neuronal 
network architecture, and the external perturbation level. 
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1. Introduction 


is 


The human brain contains about 10 11 neurons jl], and each neuron i 
connected to approximately 10 4 other neurons [ 2 J. These connections called 
synapses, are arranged in a highly complex network. They are responsible 
for neuronal communication and can be classified into two categories: elec¬ 
trical and chemical synapses [3, ij. In electrical synapses the transmission 
of information from one neuron to another is directly performed from the 
pre-synaptic cell to the post-synaptic cell via gap junctions. In chemical 
synapses, the process occurs via neurotransmitters, which cross the synaptic 
cleft and bind to receptors on the membrane of the synaptic cell [ 5 J. Neuro¬ 
transmitters may increase or decrease the probability of an action potential of 
a post-synaptic neuron, and the synapses are called excitatory or inhibitory, 
respectively [6]. Furthermore, the intensity of the chemical synapses can be 
modified, in other words, they can be minimised or potentiated [?]]. The 
mechanism responsible for these adjustments is known as synaptic plasticity 


The synaptic plasticity, that is, the ability of synapses to weaken or stren- 
then over time |9] is an important property of the mammalian brain. In 
addition, the synaptic plasticity is also related to processes of learning and 


memory 10,[llJ. This adjustment of the intensities of the chemical synapses 


can be correlated with phenomena of synchronisation of the neuronal firing 
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The occurrence of synchronisation in some specific areas of the brain may 
be associated with some diseases, such as the epilepsy and the Parkinson’s 
disease 13, 14, 0]. On the other hand, it is also responsible for some vital 
brain functions, such as processing of sensory information and motor function 

0 , [ 13 - 

Methods to suppress synchronisation have been proposed in neuroscience, 
as the introduction of external perturbations [18|, |l9|, 20]. Tass and collabo¬ 


rators have verified the possibility of desynchronisation in hippocampal neu¬ 
ronal populations through coordinated reset stimulation 22], Meanwhile, 


Popovych and collaborators have found that the introduction of a pertur¬ 
bation in a globally connected neuronal network combined with synaptic 
plasticity can provide a positive contribution to the firing synchronisation 

In this work, we study firing synchronisation in a random Hodgkin-Huxley 
neuronal network with plasticity according to spike timing-dependent plastic- 


2 

















ity (STDP). This synaptic plasticity model adjusts the connection strengths 
by means of the temporal interval between pre-synaptic and post-synaptic 


spikes [23|, 24 


Bi and Poo have reported that the change in synaptic effi¬ 
ciency after several repetitions of the experiment is due to the time difference 
of firing [25|, 26]. If one pre-synaptic spike precedes a post-synaptic spike, 


a long-term potentiation occurs, otherwise, a long-term depression appears 


A computational neuronal network specifies the connection architecture 
among neurons. A globally coupled Hodgkin-Huxley neuron model was also 
considered by Popovych and collaborators j2o[. They studied the synchro¬ 
nisation behaviour considering STDP, and found that the mean synaptic 
coupling presents a dependence on the input level. In this work, we consider 
a random neuronal network with STDP, and input, where the connections 
are associated with chemical synapses 2l[. One main result is to show that 
spike synchonisation in a neuronal network, depending on the probability of 
connections, can be improved due to spike timing-dependent plasticity. This 
improvement is also observed when an external perturbation is applied on 
the network. Another important result is that the orientation of the connec¬ 
tions among neurons with a different spike frequency affect the synchronised 
behaviour. 

This paper is organised as follows: in Section II we introduce the Hodgkin- 
Huxley neuronal model. In Section III, we show the random neuronal net¬ 
work. In Section IV, we study the synchronisation considering spike timing- 
dependent plasticity. Finally, in the last Section, we draw the conclusions. 


2. Hodgkin-Huxley neuronal network 

2.1. Hodgkin-Huxley neuronal model 

One of the most important models in computational neuroscience is the 

In this model, 


neuronal model proposed by Hodgkin and Huxley [28], [29 


the mechanism of generation of an action potential was elucidated in a series 
of experiments with the squid giant axon. They found three different ions 
currents consisting of sodium (Na), potassium (K) and leak (L) mainly due 
to chlorine. Moreover, there are voltage-dependent channels for sodium, 
potassium that control the entry and exit of these ions through the cell. The 
model is composed of a system of four coupled differential equations given 
by 

CV = I - g K n\V - E K ) - g Na m 3 h(V - E^) 
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-g L (V-E L ), 

( 1 ) 

a n (V)( 1 - n) - M v )n, 

( 2 ) 

a rn (y)( 1 - m) - f3m(V)m, 

( 3 ) 

a h (V)(l - h) - MV)h, 

( 4 ) 


where C is the membrane capacitance (measured in /jF/crn 2 ), V is the mem¬ 
brane potential (measured in mV), the function m(V) and n(V) are the 
variable of activation for sodium and potassium, and h(V) is the variable of 
inactivation for sodium. The functions a n , (3 n , a m , /3 m , (3 n are given by 
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The parameters g and E represent the conductance and reversal potentials 
for each ion. The constant / is an external current density (measured in 
pA/cm 2 ) that determines a regime of a single spike (/ = 0.0pA/cm 2 ), or 
a scheme with periodic spikes (/ = 9.0/iA/cm 2 ), as illustrated in Fig. [Qa) 
and (b), respectively. Moreover, the spikes frequency increases when the 
constant / increases. For instance, / = 9.0pA/cm 2 and / = 10.0/iA/cm 2 
approximately correspond to 67Hz and 70Hz, respectively. The parameters 
that we use in this work are presented in Table |T|! 3() . 


2.2. Network structure 

Computational models of neuronal networks depend on the architecture, 
which specifies how neurons are connected and how the dynamics is applied 
to each unit or node. In this work, we consider a random network, that 


is, the network is constructed by connecting neurons randomly |31|, [32, [33 


4 















Table 1: Parameters of the Hodgkin-Huxley neuronal model with a resting potential equal 
to — 65mV. 


Description 

Parameter 

Values 

Membrane capacity 

C 

1 yuF / cm 2 

Reversal potential for Na 

-EWa 

120 mV 

Reversal potential for K 

Ek 

-12 mV 

Reversal potential for L 

e l 

10.6 mV 

Sodium conductance 

S'Na 

120 mS/cm 2 

Potassium conductance 

9k 

36 mS/cm 2 

Leak conductance 

9l 

0.3 mS/cm 2 

External current 

I 

9.0 10.0 /iA/cni 2 



Figure 1: Dynamic firing in the Hodgkin-Huxley model, where we consider (a) / = 
0.0/rA/cm 2 that shows a single spike with a subsequent resting state, and (b) / = 
9.0/rA/cm 2 that presents a regime with periodic firing. 
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Each connection is included with probability independent from every other 
connection. Figure [2] exhibits a schematic representation of the neuronal 
network considered in this work. Each neuron is connected to others by 
randomly chosen neurons with probability p. When p = 1 we have a global 
network, where all neurons are connected. 


p<l p=l 



Figure 2: Schematic representation of a network in which each neuron is connected to 
others by randomly chosen neurons with probability p. 


We consider a random neuronal network with chemical synapses where 
the connections are unidirectional, and the local dynamics is described by 
the Hodgkin-Huxley model. The network is given by 

CVi = Ii- g K n 4 (Vi - E k ) - g Na m 3 h(Vi - £ Na ) - 

giW - E L ) + (K ~ ^ + r., ( 11 ) 

“ j -1 

where Vi is the membrane potential of neuron i (i = 1 I\ is a con¬ 

stant current density randomly distributed in the interval [9.0,10.0], to is the 
average degree connectivity, and £ij is the coupling strength from the pre- 
synaptic neuron j to the post-synaptic neuron i, that is, normally distributed 
with mean and variance equal to 0.1 and 0.02, respectively j27|. We consider 
an external perturbation Tj, so that each neuron receives an input with a 
constant intensity 7 during 1ms. This input is applied with an average time 
interval around 14ms. This value is approximately the inter-spike interval of 
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a single neuron. The neurons are excitatory coupled with a reversal potential 
V r = 20mV 20]. The post-synaptic potential s, is given by [3, 35| 


dsi _ 5(1 - Si) 


( 12 ) 


dt 1 + exp 


2.3. Spiking neurons synchronisation 

When identical neurons are coupled, the network may exhibit a complete 
synchronisation among spiking neurons, i.e., all neurons have identical time 
evolution of their action potential. We do not consider identical neurons here, 
and due to this fact, a complete synchronisation is not possible. However, a 
weak synchronisation may be observed. 

nostic of spikes synchronisation we use the order parameter given 




(13) 


where 



(14) 


where t 3 m denotes when a spike m (m = 0,1,2 ,...) of a neuron j occurs 
{tj, m < t < tj )7n+ 1). The beginning of each spike is considered when V 3 > 0. 
If the spikes times are uncorrelated, their contribution to the result of the 
summation is small [38). However, in a globally synchronised state the order 
parameter magnitude asymptotes the unity. 

Figure [3](a) shows the spiking patterns for p — 0.1, and without external 
perturbation, that is, the neuronal network presents asynchronous dynam¬ 
ics, where the points correspond to spiking neurons, and the absent points 
correspond to the resting neurons. For p = 1.0, we have globally coupled neu¬ 
rons, and we can see that the network exhibits synchronised spiking, shown 
in Fig. [3](b). The time evolution of the order parameter is plotted in Fig. 
He ). When p — 0.1 (black line), the network does not display synchronised 
behaviour, and as a result the order parameter is typically small with R fluc¬ 
tuating around 0.1. However, synchronised behaviour is observed for p = 1.0 
with order parameter values near unity (red line). 

We add an external perturbation (Tj) to analyse its effects on the syn¬ 
chronous behaviour. This way we compute the time averaged magnitude of 
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the order parameter, given by 


R 


1 


bin hi 


E 


(15) 


where the values of R have been computed by averaging over a temporal 
length of 10s after discarding a transient of 490s (hni = 490s and tfi n = 
500s). In Fig. [I] we can see the time averaged order parameter as a function 
of the probability for different input amplitudes ( 7 ). When the neuronal 
network has no input, it is possible to observe synchronised behaviour if the 
probability of connection is large enough (black circles). However, external 
inputs are able to desynchronise the spiking neurons, as shown in Fig. [4] for 
7 = 5 (red triangles) and 7 = 10 (blue squares), where the values of R are 
less than 0.9. 



Time (s) 

Figure 3: (Colour online) Space-time plots of the membrane potentials for N = 100, (a) 
p = 0.1, and (b) p = 1.0. In (c) the order parameter is calculated for p = 0.1 (black line), 
and p = 1.0 (red line). 


























Figure 4: (Colour online) Average order parameter as a function of the probability, for 
N = 100, for external perturbation with 7 = 0 (black circles), 7 = 5 (red triangles), and 
7 = 10 (blue squares). The bars represent the standard deviation. 


3. Spike timing-dependent plasticity 


One of the key principles of behavioural neuroscience is that experience 
can modify the brain structure, that is known as neuroplasticity [39j]. Al¬ 
though the idea that experience may modify the brain structure can proba¬ 


bly be traced back to the 1890s [40l . l4l|. it was Hebb who made this a central 
feature in his neuropsychological theory [42j]. 

With this in mind, we consider spike timing-dependent plasticity accord¬ 
ing to the Hebbian rule. In this plasticity the coupling strength is ad¬ 
justed based on the relative timing between the spikes of pre-synaptic and 
post-synaptic neurons 
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x f M exp(— Atij/n) , A tij >0 
13 \ ^A 2 exp(At Z j/T 2 ) , At tJ < 0 ’ 

where A Uj = t % — tj = t pos — t pie . Figure [5] exhibits the result that is obtained 
from Eq. (ITTfl) for Ai = 1.0, A 2 = 0.5, T\ = 1.8ms, and r 2 = 6.0ms. 

The initial synaptic weights £ VJ are normally distributed with mean and 
variance equal to 0.1 and 0.02, respectively. Then, they are updated accord¬ 
ing to Eq. (fl6jl . where —> e l3 + 10~ 3 Aej :/ -. In the absence of an external 
perturbation, we can verify by means of Fig. [6(a) that the averaged synaptic 
weights can be depressed (p = 0.3) or potentiated (p = 0.5, p = 0.7, and 
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Figure 5: Plasticity function (fl6l) as a function of the difference of spike timing of post- 
and pre-synaptic neuron. 



Time (s) 

Figure 6 : (Colour online) Time evolution of the averaged coupling strength for p = 0.3 
(red line), p = 0.5 (green line), p = 0.7 (blue line), and p = 1.0 (black line), where it is 
considered (a) 7 = 0, and (b) 7 = 5. 
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p — 1.0) depending on the probability of connections. If an input is applied 
on the neuronal network (Fig. [Hb), the input can have a constructive effect 
on the synaptic weights (p = 0.7 and p = 1.0) or destructive effect (p = 0.3 
and p = 0.5), depending on the probability of connections. 



Figure 7: (Colour online) Average order parameter versus probability of connection for 
the cases with (black circles) and without (red triangles) STDP, where we consider (a) 
7 = 0, (b) 7 = 3, (c) 7 = 5, and (d) 7 = 10. The bars represent the standard deviation. 

The time averaged order parameter in terms of the probability of connec¬ 
tion is showed in Fig. [71 considering with (black circles) and without (red 
triangles) STDP. I 11 the case for without external perturbation (Fig. [7^) we 
can see that the values of R without STDP are less than with STDP, namely 
STDP is producing a positive effect on the synchronisation. Increasing the 
amplitude of the external perturbation, without spike timing-dependent plas¬ 
ticity, the desynchronisation is induced in the neuronal network (Fig. [I]). 
However, considering STDP in the perturbed network is possible to observe 
alterations to the dynamic behaviour in relation to synchronised states. Fig¬ 
ure [71(b) shows that the STDP enhances the synchronisation for p approxi¬ 
mately greater than 0.5 due to a constructive effect on the dynamics of the 
synaptic weights. On the other hand, for p less than 0.5 the STDP decreases 
the values of the time averaged order parameter, as a result of depressed 
synaptic weights. Increasing the input intensity (7 = 5), it is possible to 
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verify synchronisation when the network has neuroplasticity (Fig. [7fc). We 
can also see an abrupt transition from a desynchronised to a synchronised 
regime. When the input intensity increases, it is necessary to increase the 
probability of connection for the STDP counteracts the suppression of syn¬ 
chronisation. Figure[7](d) exhibits a situation such that the synchronisation is 
only obtained when the neuronal network presents a global coupling (p — 1). 
With a strong input, the STDP does not lead the network to a potentiation 
of synaptic weights, and this way the synchronisation is suppressed by an 
external input. 




Figure 8 : (Colour online) (a) Absolute value of the plasticity function (fl 6 l) as a function 
of the difference of spike timing of post- and pre-synaptic neuron, where the curves are 
the potentiation (red line) and the depression (black line), (b) Average time difference 
versus probability of connections for 7 = 0 (red circles), 7 = 3 (blue squares), 7 = 5 (black 
triangles), and 7 = 10 (green diamonds), (c) Critical probability p c as a function of the 
input level 7 . The linear fit is given by the equation p c = O.O 67 + 0.37. 


The spike synchronisation depends on the probability of connections p in 
a way showing an abrupt transition. There is a critical point for p = p c that 
can be found by means of the intersection between the curves of potentiation 
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and depression. Figure |SJa) exhibits the point of intersection with value of 
A tij approximately equal to 1.8. For > 1.8 the depression (red line) 
of the synaptic strength is larger than the potentiation (black line), while 
that for Atij <1.8 the potentiation is larger than the depression. With this 
value of At^ we can obtain p c plotting A tij as a function of p, as shown in 
Fig. [51(b) . It can be seen that the value of At t] decays, and when cross the 
value Atij ~ 1-8 we have the critical value of the probability of connections. 
Moreover, when p increases not only Adecreases, but also the variance of 
the inter-spike intervals decreases. Then, we compute the critical probability 
as a function of the input level (Fig. IHt), where we verify a linear increase 
given by the equation p c = O.O 67 + 0.37. 

Figure EKc) shows a discontinuous transition between the synchronised 
and the desynchronised regime. Our aim is to understand how this discon¬ 
tinuous transition appears when the probability of connections is varied. For 
this reason, we build a network according to a schematic representation that 
is showed in Fig. |9j The scheme represents a network of neurons with high 
(cyan ball) and low (yellow ball) spiking frequency. The red arrows represent 
the connections between neurons from high to low frequency, while the black 
arrows represent the connections from low to high frequency. 



Figure 9: (Colour online) Schematic representation of a network of neurons with high 
(cyan ball) and low (yellow ball) spiking frequency. The red arrows represent the connec¬ 
tions between neurons from high to low frequency, while the black arrows represent the 
connections from low to high frequency. 


Based on the schematic representation that is illustrated in Fig. [9l we 
consider a neuronal network with N = 100, p = 0.47, and 7 = 3 . We separate 
the neuronal network into 50 neurons with high frequency (values of /; within 
the range [9.0,9.!]) and 50 neurons with low frequency (values of Jj within 
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the range [9.9,10]). Figure [TU] exhibits the time averaged coupling strength as 
a function of the percentage of connections from neurons with high frequency 
to neurons with low frequency. We can see that the time average coupling 
strength depends on the connections. Considering the case for a small per¬ 
centage of connections from HFN to LFN, the average coupling is small, 
indicating absence of synchronisation, a situation that changes with increas¬ 
ing the percentage of connections to a synchronised state. This means that, 
when the coupling strengths increase, a desynchronised state can suddenly 
become synchronised. Consequently, the abrupt transition from desynchro¬ 
nised to synchronised state, that is observed in Fig. [7] is due to directed 
synapses among spiking neurons with high and low frequency. 



Figure 10: Time averaged coupling strength versus percentage of connections from neurons 
with high frequency (HFN) to neurons with low frequency (LFN). We consider p = 0.47 
and 7 = 3. 


4. Conclusion 

We have been studying a neuronal network model with spiking neurons. 
We have chosen, as local dynamics, the Hodgkin-Huxley model due to the 
fact that it has essential features of spiking dynamics. The Hodgkin-Huxley 
model is a coupled set of nonlinear differential equations that describes the 
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ionic basis of the action potential. These equations are able to reproduce 
biophysical properties of the action potential. 

We have used a random coupling architecture where the connections are 
randomly distributed according to a probability. When the probability is 
equal to unity we have a globally coupled network. The connections were 
considered unidirectional representing excitatory chemical synapses. 

We have studied the effects of spike timing-dependent plasticity on the 
synchronisation in a Hodgkin-Huxley neuronal network. Studies about spike 
synchronisation are important to understand not only progressively degener¬ 
ative neurological disorders, but also processing of the sensory information. 
Popovych and collaborators 20] showed that STDP combined with an exter¬ 
nal perturbation can improve the spike synchronisation in a globally neuronal 
network. The novelty in this paper is that we have considered a random neu¬ 
ronal network and we have verified that the spike synchronisation depend 
on the probability of connections. Considering a strong external perturba¬ 
tion the spike synchronisation is suppressed. However, when there is STDP, 
depending on the probability of connections, the synchronisation in the per¬ 
turbed network can be improved due to a constructive effect on the synaptic 
weights. 

We have also shown that the direction of synapses has an important role 
on the effects of spike timing-dependent plasticity on the synchronisation in 
a random Hodgkin-Huxley neuronal network. 
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